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FIRST ORDER KINETICS 


DECUS Program Library Write-up 


DECUS No. 8-133 


This program is designed to calculate the rate constants and 


error parameters for first order rate processes. The program 
provides options for plotting the deviation between observed and 
calculated quantities on a CRT and for iterating to a "best" infinity 
value. The program may also be used in other cases in which one 
wishes to correlate the In of one quantity with another, as in linear 
free energy relationships. 

Specifications: The program uses most of memory for the main 
program, the subroutines and data storage. A storage map is: 



63-1705 J Main P ro g ra ™ 


2000-2177 PLOT Subroutine (attached) 

3400-3666 LSQ Subroutine (DECUS 8 

4000-4057 MESAGE Subroutine (DIGITAL-8-18-U mod. 


(see below) 


4075-4122 DECPRT Subroutine (DIGITAL-8-22-U) 
4200-4311 SICONV Subroutine (DIGITAL-8-28-U) 
4315-4366 FIX and FLOAT Subroutines (Cf Floating Point 


System Manual) 


4400-4534 MINMAX Subroutine (attached) 

4757-7577 Floating Point Interpreter "C" (8-5C-S) 

2200-2377 ABS array (holds absorbance values) 
2400-2577 DIF array (holds devations between observed 


and calculated y's) 


2600-2777 CALC array (holds calculated values of y's) 
3000-3177 Y array (holds observed y values) 

3200-3377 T array (holds values of time) 


The MESAGE subroutine was modified by replacing the 


TLS 
TSF 
JMP. -1 


printing sequence with 


TSF 
JMP. -1 
TLS 
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which is compatible with the typeout routine in the floating point 
interpreter. 

The program is written for a 4K PD P-8. If an oscilloscope 
display (Type 34D) is available, it may be used to examine the 
deviations between observed and calculated quantities. If it is not 
available, the program will still operate normally. 

Description: 

The program is loaded in the normal fashion using the binary 
loader. If the binary tape which is supplied is used, it will stop 
half way through (after loading the floating point interpreter). The 
CONTINUE key is then depressed in order to load the remainder of 
the program. The binary tape contains the main program and all 
needed subroutines. 

The starting address is 200. When started, the computer 
will type: 


FIRST ORDER KINETICS 
DATA IN ABSORBANCE, M=1 
DATA IN TRANSMITTANCE, M=2 
DATA IN VOLUME, M=3 
M= 

The operator then types the proper number and return. Then appears 

N IS THE NUMBER OF DATA POINTS 
N= 

The operator types the number of points. Then: 

THE INFINITY VALUE = 

The operator types the value of the measured quantity which is expected 
at infinite time. Then follows 


NOW TYPE observed values AND TIMES 

The operator types the observed quantity and the corresponding time, 

and does this N times. As a check, when the last point has been accepted, 
there is typed: * 


ALL data have been ACCEPTED 
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If M=1 or 3, the rate constant is calculated by the subroutine 
LSQ (cf. DECUS 8 ) using 

yj_ = In (ABSj-ABSoj) = mt£ + b 

The rate constant and other quantities are typed as well as the calculated 
and observed y's and the deviation, (cf. the following output sheet. Fig 1) 
If M=2, there is typed 

100 PC transmittance= 

The operator types the value of the measured quantity when no light 
absorption occurs. Then follows: 

O PC TRANSM1TTANCE= 


The dark current is then typed. If the former quantity is called MAX 
and the latter MIN, the absorbance values are given by 


ABSi= 


1 

2. 303 


LN 


( MAX-MIN \ 
TOBS-MIN/ 


where TOBS is the observed value of transmittance at the given time. 
The remainder of the calculation proceeds as before. 


After the data have been typed, there appears 
IF PLOT IS DESIRED, RAISE KEY 1, PRESS CONTINUE 


If the oscilloscope display is attached and one wished to examine the 
deviation as a function of time, switch register key 1 is raised and 
CONTINUE is depressed. The display will continue as long as key 1 
remains up. 

The plot is formed by first searching for the maximum and 
minimum values in the arrays using the subroutine MINMAX which is 
a modification of our standard version. Having the minimum and 
maximum values, the subroutine PLOT scales all values to be in the 
range 0-1777 and then plots the x and y axes, and the points in the DIF 
array. This is a particularly valuable option for it permits the user 
to quickly see if his data fit a first order process. One which fits 
well gives the type of scatter seen in Fig 2a, whereas one which does 
not (because it is not a first order reaction or because a poor infinity 
value was chosen) gives a curve of the type seen in Fig 2b. 
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When key 1 is dropped, there appears 


ITER IS THE NUMBER OF ITERATIONS, ITER= 

If the value -1 is given, the computer halts. Pressing CONTINUE restarts 
the program at the beginning. If the value 0 is given, the computer will 
print out the observed quantities (absorbance or volume) and the values 
obtained from the calculated y's. The computer then halts and may be re¬ 
started by depressing the CONTINUE key. If a positive integer is given, 
there will be typed: 

MULTIPLY VINF BY 

The operator types a value by which the infinity value is to be multiplied in 
the first iteration step. Then appears 

VINF= 


The operator types the infinity value to be used (usually the one given previously. 
If MODE=2, the infinity value must be given in absorbance rather than trans¬ 
mittance. The infinity absorbance may be obtained by first using ITER=0. After 
the data have been typed, depressing CONTINUE will again permit a plot to be 
observed and iterations to be carried out. 

The computer will perform the desired number of iterations. Whenever 
the rms deviation is found to increase, the fractional change in the multiplier 
is halved, and the reciprocal is taken for the next iteration step. Thus, if the 
value typed above were 1. 1 (giving a 10% increase in VINF in the first iteration), 
and the rms deviation were found to increase, the next step would use as a 
multiplier 1/ 1. 05 or 0. 952. This multiplier would continue to be used until the 
rms deviation no longer decreases. Then, the multiplier is again changed and 
will be 1. 024. 

This is not the ideal way in which to iterate to the "best" infinity value 
since the least squares criterion remains the minimization of the error in y, 
not in the observed quantity (absorbance or volume). However, a better cal¬ 
culation is not possible with the storage available in a 4K computer, and in our 
experience, the difference in result between the present calculation and the 
"correct" method is normally negligible. 

The program permits the use of up to 42 data points. A larger number 
of points will lead to malfunction. If a double space is desired between lines 
of output data, the following locations may be changed: 
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Location 

Now 

Change to 

55 

0000 

7777 

1565 

7200 (CLA) 

7240 (CLA CMA) 

1244 

7200 (CLA) 

7240 (CLA CMA) 

1236 

7200 (CLA) 

7240 (CLA CMA) 
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FIRST ORDER KINETICS 
DATA IN ABSORBANCE/ M = 1 
DATA IN TRANSMITTANCE/ M = 2 
DATA IN VOLUME/ M = 3 
M = J_ 

N IS THE NUMBER OF DATA POINTS 
N = L2 

THE INFINITY VALUE =0.135 

NOW TYPE OBSERVED VALUES AND TIMES 

0.700 0 

0.665 To 

0-626 2£ 

0.598 30 
0.568 40 
0.510 60 
0.463 80 
0.425 100 
0.373 140 
0.304 180 
0.245 240 
0.200 320 

ALL DATA HAVE BEEN ACCEPTED 
THE RATE CONSTANT IS + 0.67 42559 E-02 
THE INTERCEPT IS -0.5654722E+00 
THE RMS ERROR IS +0.2329548E-0 1 
R = +0.999356 5E+00 

THE PERCENT ERROR IN K IS +0.6703959E+00 

DIF 

+0.5457162E-02 
+ 0. 1979947E-02 
+0.1098728E-01 
+ 0.2278 68 5E-02 
+ 0. 18 42617E-02 
+ 0. 1080358 E-0 1 
+ 0.986409 IE-02 
-0. 18539 42E-02 
“0.7394647E-01 
-0.1276254E-02 
+ 0.23588 18E-01 
+ 0. 1027 536E-01 


T 

+ 0.0000000E+00 
+ 0. 1000000E+02 
+ 0.2000000E+02 
+ 0.3000000E+02 
+0.4000000E+02 
+0.6000000E+02 
+0.8000000E+02 
+0.1000000E+03 
+0.1400000E+03 
+0.1800000E+03 
+0•2400000E+03 
+0•3200000E+03 


Y 

-0. 5709294E+00 
- 0•6 348778 E+00 
-0.7113107E+00 
-0.7700277E+00 
-0.8370172E+00 
-0.9808293E+00 
-0.1114741E+01 
-0. 123787 4E+01 
-0.1435484E+01 
-0.1777856E+01 
-0.2207275E+01 
-0.2733367E+01 


Y CALC 

-0.56 54722E+ 0 0 
-0.6328978E+00 
-0.7003234E+00 
-0.7677490E+00 
-0.8351746E+00 
-0.9700257E+00 
-0. 1 104876E+01 
-0. 1239728E+0 1 
-0. 1 509430E+0 1 
-0.1779132E+01 
-0.2 18 3686E+0 1 
-0.2723091E+01 


IF PLOT IS DESIRED/ RAISE KEY 1/ PRESS CONTINUE 


Fig. la Typical run using mode 1 (data in absorbance). The program 

continues after that shown by asking for the number of iterations 
The quantities typed by the operator are underlined. 
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FIRST ORDER KINETICS 
DATA IN ABSORBANCE, M = 1 

DATA IN TRANSMITTANCE, M = 2 
DATA IN VOLUME, M = 3 
M = 2_ 

N IS THE NUMBER OF DATA POINTS 
N = _6_ 

THE INFINITY VALUE =0.8328 

NOW TYPE OBSERVED VALUES AND TIMES 

.2995 0. 

.3366 20 
.3704 40 
.4090 60 
.47 58 100 

.7310 320 


ALL DATA HAVE BEEN ACCEPTED 
100 PC TRANSMITTANCE = KJ 
0 PC TRANSMITTANCE = 0j_l 
THE RATE CONSTANT IS +0.6754450E- 
THE INTERCEPT IS -0.5707473E+00 
THE RMS ERROR IS +0.48 7 1 432E-02 
R = +0.9999778E+00 

THE PERCENT ERROR IN K IS +0.1351 


+0.0000000E+00 
+0.2000000E+02 
+0.4000000E+02 
+0.6000000E+02 
+0.1000000E+03 
+0.3200000E+03 


-0.5708 574E+00 
-0.7 1 13739 E+00 
-0.8370701E+00 
-0.9807610E+00 
-0. 1237778E+01 
-0.2 7 3 4 0 48 E+0 1 


02 


247 E+00 
Y CALC 

-0.5707473E+00 
- 0.7058 364E+00 
-0.8409254E+00 
-0.9760143E+00 
-0.1246192E+01 
-0.27 32 171 E+01 


DIF 

+0.1100301E-03 
+ 0.5537 509 E-02 
-0.38 55347E-02 
+ 0.47 4667 5E-02 
-0.841379 IE-02 
+ 0. 18768 31 E-02 


IF PLOT IS DESIRED, RAISE KEY 1, PRESS CONTINUE 


ITER IS THE NUMBER OF ITERATIONS, ITER = _0 
THE RATE CONSTANT IS +0.6754450E-02 
INTERCEPT IS -0.5707473E+00 
RMS ERROR IS +0.487 1432E-02 
+0.9999778E+00 

PERCENT ERROR IN K IS +0.1351247E+00 

RMS ERROR BASED ON OBSERVED QUANT. IS +0.1794132E-02 
INFINITY VALUE IS +0.1 350 1 41 E+00 


THE 
THE 
R = 

THE 
THE 
THE 

TIME 

+ 0.0000000 E+ 0 0 
+0•2000000E+02 
+0.4000000E+02 
+ 0.6000000E+02 
+ 0.1000000E+03 
+0.3200000E+03 


ABSORBANCE 
+0.7000547E+00 
+ 0.62598 32E+00 
+0.5679912E+00 
+0.5100395E+00 
+0.4250417E+00 
+ 0. 1999698 E+00 


CALC 

+0.7001166E+00 
+0.6287096E+00 
+0.5663254E+00 
+ 0.51 18238 E+00 
+0.4226120E+00 
+ 0.2000919 E+00 


DIF 

-0.618 69 62 E-04 
-0.2726435E-02 
+ 0. 16658 30E-02 
-0. 178 4324E-02 
+0.2429723E-02 
-0.1220405E-03 


Fig. lb Typical run using mode 2 (data in transmittance). 
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Fig. 2. Computer drawn plots of difference between observed 
and calculated Quantities for a first order reaction. The 
upper plot corresponds to the data in Fig. l and shows a 
normal scatter. The lower plot is for data which are not 
well represented as a first order process. 
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/FIRST ORDER KINETICS 

/PAGE 0 

CALL= 4400 

FLPT= 4407 

RETURN=5400 

FSQRT=2 

FEXP=6 

FLOG= 7 

INPUT=13 

OUTPUTS 14 




*7 



0007 

5600 


5600 


0010 

0000 


0 


001 1 

0000 

IT, 

0 

/INDEX REGISTERS 

0012 

0000 

IY, 

0 


0013 

0000 

I CALC, 

0 


0014 

0000 

IDIF, 

0 


0015 

0000 

I ABS, 

*55 

0 


0055 

0000 

*67 

0 


0067 

0000 

CORR, 

FLTG 0.0 


0070 

0000 




0071 

0000 


/CORRECTION 

FACTOR 

0072 

0000 

TEST, 

FLTG 0.0 


0073 

0000 




0074 

0000 




0075 

0000 

FACTOR, 

FLTG 0.0 


0076 

0000 




0077 

0000 




0 100 

0000 

S, 

FLTG 0.0 


0 101 

0000 




0 102 

0000 


/SLOPE 


0 103 

0000 

B, 

FLTG 0.0 


0104 

0000 




0 105 

0000 


/INTERCEPT 


0 106 

0000 

R, 

0 

/CORRELATION COEFF 

0 107 

0000 


0 


0 1 10 

0000 


0 


0111 

0000 

ERROR, 

FLTG 0.0 


0 1 12 

0000 




0 113 

0000 


/ERROR 


0 114 

0000 

G, 

FLTG 0.0 


0115 

0000 




0 116 

0000 


/FLOAT N 
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0 1 17 

0000 

N, 

0 

/NUMBER OF 

0 120 

3200 

LT, 

3200 

/LOCC T) 

0 121 

3000 

LY , 

3000 

/LOCCY) 

0 122 

2600 

LCALC, 

2600 

/LOCC CALC) 

0 123 

2400 

LDIF, 

2400 

/LOCC DIF) 

0 124 

2200 

LABS, 

2200 

/LOCCABS) 

0 125 

0000 

M, 

0 

/MODE 

0 126 

0000 

FIN, 

FLTG 

0.0 

0 127 

0000 




0 130 

0000 







/FINAL VALUE 

0131 

4000 

MESAGE, 

4000 


0 132 

4075 

DECPRT, 

4075 


0 133 

4200 

SICONV, 

4200 


0 134 

4315 

FIX, 

4315 


0 135 

4355 

FLOAT, 

4355 


0 136 

3400 

LSQ, 

3400 


0 137 

2000 

PLOT, 

2000 


0140 

4400 

MINMAX, 

4400 


0 141 

0000 

ITEk, 

0 


0142 

0000 

NEGN, 

0 


0 143 

7776 

M2, 

-2 


0 144 

7777 

CONST, 

FLTG 

0.434294 

0145 

3362 




0 146 

6724 




0 147 

0001 

ONE, 

1 


0150 

0000 

ZERO, 

FLTG 

0.0 

0151 

0000 




0152 

0000 




0153 

0000 

FLAG, 

0 


0154 

0000 

XLOC, 

0 

/CONSTANTS 

0155 

0000 

YLOC, 

0 


0156 

0000 

XMIN, 

FLTG 

0.0 

0 157 

0000 




0160 

0000 




0161 

0000 

XRNG, 

FLTG 

0.0 

0162 

0000 




0163 

0000 




0164 

0000 

YMIN, 

FLTG 

0.0 

0165 

0000 




0166 

0000 




0167 

0000 

YRNG, 

FLTG 

0.0 

0 170 

0000 




0171 

0000 




0172 

0000 

XZERO, 

FLTG 

0.0 

0 173 

0000 




0 174 

0000 




0175 

0000 

YZERO, 

0 


0176 

0000 


0 


0177 

0000 


0 



POINTS 
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0200 

7200 

/PAGE 1 
*200 
START, 

0201 

3141 


0202 

6032 


0203 

6046 


0204 

4531 


0205 

4543 

TEXT /%# 

0206 

061 1 

FI 

0207 

2223 

RS 

0210 

2440 

T 

0211 

1722 

OR 

0212 

0405 

DE 

0213 

2240 

R 

0214 

1311 

K I 

0215 

1605 

NE 

0216 

241 1 

TI 

0217 

0323 

CS 

0220 

4543 

%# 

0221 

0401 

DA 

0222 

2401 

TA 

0223 

4011 

I 

0224 

1640 

N 

0225 

0102 

AB 

0226 

2317 

SO 

0227 

2202 

RB 

0230 

01 16 

AN 

0231 

0305 

CE 

0232 

5440 

, 

0233 

1540 

M 

0234 

7540 

= 

0235 

6100 

1/ 

0236 

4531 


0237 

4543 

TEXT /%# 

0240 

0401 

DA 

0241 

2401 

TA 

0242 

401 1 

I 

0243 

1640 

N 

0244 

2422 

TR 

0245 

01 16 

AN 

0246 

2315 

SM 

0247 

1 124 

I T 

0250 

240 1 

TA 

0251 

1603 

NC 

0252 

0554 

E, 

0253 

4015 

M 

0254 

407 5 

= 

0255 

4062 

2 

0256 

4543 



OF FIRST ORDER KINETICS PROGRAM 
CL A 

DCA ITER /SET FOR NO ITERATIONS 

KCC 

TLS /PREPARE FOR OUTPUT 

CALL MESAGE 


CALL MESAGE 
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0257 

040 1 

DA 


0260 

240 1 

TA 


0261 

401 1 

I 


0262 

1640 

N 


0263 

2617 

VO 


0264 

1425 

LU 


0265 

1505 

ME 


0266 

5440 

» 


0267 

1540 

M 


0270 

7540 

= 


027 1 

6300 

3/ 


0272 

4531 


CALL MESAGE 

0273 

4543 

TEXT /%# 


0274 

1540 

M 


0275 

7 540 

= 


0276 

0000 

/ 


0277 

4533 


CALL SICONV 

0300 

3125 


DCA M 

0 301 

4531 


CALL MESAGE 

0 302 

4543 

TEXT /%# 


0303 

1640 

N 


0304 

1 123 

I S 


0 30 5 

4024 

T 


0306 

1005 

HE 


0 307 

4016 

N 


0310 

2515 

UM 


031 1 

0205 

BE 


0312 

2240 

R 


0313 

1706 

OF 


0314 

4004 

D 


0315 

0124 

AT 


0316 

0140 

A 


0317 

2017 

PO 


0 320 

1 116 

IN 


0 321 

2423 

TS 


0 322 

4543 

%# 


0 323 

1640 

N 


0 324 

7 540 

= 


0325 

0000 

/ 


0 326 

4533 


CALL SICONV 

0 327 

3117 


DCA N 

0330 

4531 


CALL MESAGE 

0331 

4543 

TEXT /%# 


0332 

2410 

TH 


0333 

0540 

E 


0334 

1 1 16 

IN 


0335 

061 1 

FI 


0336 

1611 

NI 


0 337 

2431 

TY 
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0 340 

4026 

V 

0341 

01 14 

AL 

0 342 

2 50 5 

UE 

0343 

407 5 

= 

0344 

0000 

/ 

0345 

4407 


0346 

0013 


0 347 

6 126 


0350 

0000 


0351 

4531 


0352 

4543 

TEXT 

0353 

1617 

NO 

0354 

2740 

W 

0355 

2431 

TY 

0356 

2005 

PE 

0357 

4017 

0 

0360 

0223 

BS 

0361 

0 522 

ER 

0362 

2605 

VE 

0363 

0440 

D 

0364 

2601 

VA 

0365 

1425 

LU 

0366 

0523 

ES 

0367 

4001 

A 

0370 

1604 

ND 

0 37 1 

4024 

T 

0372 

1115 

IM 

0373 

0523 

ES 

0374 

4543 

%# 

0375 

0000 

/ 

0376 

5777 


0377 

0400 

PG2, 


FLPT 
INPUT 
FPUT FIN 
FEXT 

CALL MESAGE 

/%# 


JMP I PG2 
400 
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/PAGE 2 OF FIRST ORDER 


0400 

7200 

*400 

MIN=167 

MAX= 172 
TEMP= 175 

CL A 

0401 

1 1 17 


TAD N 

0402 

7041 


CIA 

0403 

3142 


DCA NEGN 

0404 

1 124 


TAD LABS 

0405 

3015 


DCA IABS 

0406 

1 120 


TAD LT 

0407 

3011 


DCA IT 

0410 

4407 

READ, 

FLPT 

0411 

0013 


INPUT 

0412 

6415 


FPUT I I ABS 

0413 

0013 


INPUT 

0414 

641 1 


FPUT I IT 

0415 

0000 


FEXT 

0416 

2142 


ISZ NEGN 

0417 

5221 


JMP .+2 

0420 

5230 


JMP NEXT 

0421 

2015 


ISZ IABS 

0422 

2015 


ISZ IABS 

0423 

2015 


ISZ I ABS 

0424 

201 1 


ISZ IT 

0 425 

201 1 


ISZ IT 

0 426 

2011 


ISZ IT 

0427 

5210 


JMP READ 

0 430 

4531 

NEXT, 

CALL MESAGE 

0431 

4543 

TEXT /%# 


0432 

0114 

AL 


0433 

1440 

L 


0434 

0401 

DA 


0435 

2401 

TA 


0436 

4010 

H 


0437 

0126 

A V 


0440 

0540 

E 


0441 

0205 

BE 


0442 

0516 

EN 


0443 

4001 

A 


0444 

0303 

CC 


0445 

0520 

EP 


0446 

2405 

TE 


0 447 

0400 

D/ 


0450 

7200 


CLA 

0451 

1 125 


TAD M 

0452 

1143 


TAD M2 

0453 

7440 


SZA 

0454 

5363 


JMP YCALC 

0455 

7200 


CLA 

0456 

1 124 


TAD LABS 


KINETICS PROGRAM 


/COMPLEMENT N 


/INITIALIZE REGISTERS 


/READ ABS 
/AND STORE 
/READ TIME 
/AND STORE 

/INCREMENT COUNTER 
/CONTINUE LOOP 
/CONTINUE PROGRAM 


/LOOP 


/TEST MODE 

/IS M = 2? 

/NO: CALC Y'S 

/YES: CONV. TRANS TO ABS 


14 





0457 

3015 


DCA IABS 


0460 

4531 


CALL MESAGE 


046 1 

4543 

TEXT /%# 



0462 

6160 

10 



0463 

6040 

0 



0464 

2003 

PC 



0465 

4024 

T 



0466 

2201 

RA 



0467 

1623 

NS 



0470 

151 1 

MI 



047 1 

2424 

TT 



0472 

01 16 

AN 



0 47 3 

0305 

CE 



0 47 4 

4075 

= 



0 47 5 

4000 

/ 



0476 

4407 


FLPT 


0477 

0013 


INPUT 

/ 100% TRANS 

0 500 

6172 


FPUT MAX 


0 501 

0000 


FEXT 


0 502 

4531 


CALL MESAGE 


0503 

4543 

TEXT /%0 



0 504 

6040 

0 



0 505 

2003 

PC 



0 506 

4024 

T 



0 507 

2201 

RA 



0510 

1623 

NS 



051 1 

151 1 

MI 



0512 

2424 

TT 



0513 

01 16 

AN 



0514 

0305 

CE 



0515 

4075 

= 



0516 

4000 

/ 



0517 

4407 


FLPT 


0 520 

0013 


INPUT 

/DARK CURRENT 

0 521 

6167 


FPUT MIN 


0 522 

5172 


FGET MAX 


0 523 

2167 


FSUB MIN 


0 524 

6172 


FPUT MAX 


0 525 

0000 


FEXT 


0 526 

7200 


CL A 


0 527 

1117 


TAD N 


0 530 

7041 


CIA 


0 531 

3142 


DCA NEGN 


0 532 

4407 

LOOP1 » 

FLPT 


0533 

5415 


FGET I IABS 

/TRANSMI TTANCE 

0534 

2167 


FSUB MIN 

/-DARK CURRENT 

0535 

6175 


FPUT TEMP 


0 536 

5172 


FGET MAX 


0 537 

4175 


FDI V TEMP 
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0540 

0007 


FLOG 


0 541 

3144 


FMPY 

CONST 

0 542 

6415 


FPUT 

I I ABS 

0543 

0000 


FEXT 


0544 

2015 


ISZ 

I ABS 

0 545 

2015 


I SZ 

I ABS 

0 546 

2015 


I SZ 

I ABS 

0 547 

2142 


I SZ 

NEGN 

0550 

5332 


JMP 

LOOP 1 

0551 

4407 


FLPT 


0 552 

5126 


FGET 

FIN 

0553 

2167 


F SUB 

MIN 

0 554 

6126 


FPUT 

FIN 

0555 

5172 


FGET 

MAX 

0556 

4126 


FDI V 

FIN 

0557 

0007 


FLOG 


0 560 

3144 


FMPY 

CONST 

0 561 

6126 


FPUT 

FIN 

0 562 

0000 


FEXT 


0563 

7200 

Y CALC# 

CL A 


0564 

1 124 


TAD 

LABS 

0565 

3015 


DCA 

I ABS 

0 566 

1 121 


TAD 

LY 

0567 

3012 


DCA 

IY 

0570 

1 120 


TAD 

LT 

0 571 

301 1 


DCA 

IT 

0572 

1 122 


TAD 

LCALC 

0 57 3 

3013 


UCa 

I CALC 

0574 

1 123 


TAD 

LDIF 

0 57 5 

3014 


DCA 

IDIF 

0576 

5777 


JMP 

I PG3 

0577 

0600 

PG3, 

600 



/LN (T’/T) 

/LOG = ABS 
/STOKE hBS 

/CONTINUE LOOP 
/TRANSFORM INFINITY VALUE 


/SET UP REGISTERS 
/FOR CALC OF Y 


<9 
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/PAGE 3 OF FIRST ORDER 

0600 

7200 

F ACTOR= 7 5 
TEM2= 1 56 
*600 

CL A 

0601 

1 1 17 


TAD N 

0 602 

7041 


CIA 

0603 

3142 


DCA NEGN 

0604 

4407 


FLPT 

0605 

5150 


FGET ZERO 

0606 

6111 


FPUT ERROR 

0607 

6114 


FPUT G 

06 10 

5415 


FGET I IABS 

0611 

2126 


FSUB FIN 

0612 

6156 


FPUT TEM2 

0613 

0000 


FEXT 

0614 

7200 


CL A 

0615 

1 157 


TAD TEM2+1 

0616 

7510 


SPA 

0617 

5227 


JMP NEG 

0620 

7200 


CLA 

0621 

1 147 


TAD ONE 

0622 

4535 

A1 > 

CALL FLOAT 

0623 

4407 


FLPT 

0624 

6075 


FPUT FACTOR 

0625 

0000 


FEXT 

0626 

5233 


JMP L00P3 

0627 

7200 

NEG* 

CLA 

0630 

1 147 


TAD ONE 

0631 

7041 


CIA 

0632 

5222 


JMP A1 

0633 

4407 

L00P3* 

FLPT 

0634 

5415 


FGET I I ABS 

0635 

5415 


FGET I I ABS 

0636 

2126 


FSUB FIN 

0637 

3075 


FMPY FACTOR 

0640 

0007 


FLOG 

0641 

6412 


FPUT I IY 

0 642 

0000 


FEXT 

0643 

2012 


ISZ IY 

0644 

2012 


ISZ IY 

0645 

2012 


ISZ IY 

0646 

2015 


ISZ IABS 

0647 

2015 


ISZ IABS 

0 6 50 

2015 


ISZ IABS 

0651 

2142 


ISZ NEGN 

0 6 52 

5233 


JMP L00P3 

0653 

4536 


CALL LSD 

0654 

4531 


CALL MESAGE 

0655 

4543 

TEXT /%# 


0656 

2410 

TH 



KINETICS PROGRAM 


/DO VALUES INCREASE 
/OR DECREASE? 


/SET FACTOR APPROPRIATELY 


/CALCLATE Y 
/STORE Y 
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0657 

0540 

E 

G 660 

2201 

RA 

066 1 

2405 

TE 

0662 

4003 

C 

0663 

17 16 

ON 

0 664 

2324 

ST 

0665 

0116 

AN 

0666 

2440 

T 

0667 

1 123 

I S 

0670 

4000 

/ 

0671 

4407 


0672 

5150 


0673 

2100 


0674 

0014 


0675 

0000 


0676 

4531 


0677 

4543 

TEXT 

0700 

2410 

TH 

0701 

0540 

E 

0702 

1 116 

IN 

0703 

2405 

TE 

0704 

2203 

RC 

0705 

0520 

EP 

0706 

2440 

T 

0707 

1 123 

IS 

0710 

4000 

/ 

0711 

4407 


0712 

5103 


07 13 

0014 


07 14 

0000 


0715 

4531 


07 16 

4543 

TEXT 

0717 

2410 

TH 

0720 

• 0540 

E 

0721 

2215 

RM 

0722 

2340 

S 

0723 

0522 

ER 

0724 

2217 

RO 

0725 

2240 

R 

0726 

1 123 

IS 

0727 

4000 

/ 

0730 

4407 


0731 

511 1 


0 7 32 

0014 


0733 

0000 


0734 

4531 


0735 

4543 

TEXT 

0736 

2240 

R 

0737 

7 540 

2 


FLPT 

FGET ZERO 
F SUB S 
OUTPUT 
FEXT 

CALL MESAGE 

/%£ 


FLPT 
FGET B 
OUTPUT 
FEXT 

CALL MESAGE 

/%# 


FLPT 

FGET ERROR 

OUTPUT 

FEXT 

CALL MESAGE 

/%# 
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0740 

0000 

/ 

0741 

4407 


0 7 42 

5106 


0743 

0014 


0744 

0000 


0745 

4531 


0746 

4543 

TEXT /%# 

0747 

2410 

TH 

0750 

0540 

E 

0751 

2005 

PE 

07 52 

2203 

RC 

0 7 53 

0516 

EN 

0754 

2440 

T 

0755 

0522 

ER 

0 7 56 

2217 

RO 

0757 

2240 

R 

0760 

1 1 16 

IN 

0761 

4013 

K 

0762 

401 1 

I 

0763 

2340 

S 

0764 

0000 

/ 

0765 

4407 


0766 

5156 


0767 

0014 


0770 

0000 


0771 

5772 


0772 

1000 

PG4, 




FLPT 
FGET R 
OUTPUT 
FEXT 

CALL MESAGE 


FLPT 

FGET TEM2 

OUTPUT 

FEXT 

JMP I PG4 
1000 
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/PAGE 4 OF FIRST ORDER 
*1000 


1000 

7200 


CL A 


1001 

1 141 


TAD 

ITER 

1002 

7440 


SZA 


1003 

5777 


JMP 

ITRATE+1 

1004 

7200 


CL A 


1005 

1 117 


TAD 

N 

1006 

7041 


CIA 


1007 

3142 


DCA 

NEGN 

1010 

1 120 


TAD 

LT 

1011 

301 1 


DCA 

I T 

10 12 

1 121 


TAD 

LY 

1013 

3012 


DCA 

IY 

1014 

1 122 


TAD 

LCALC 

1015 

3013 


DCA 

I CALC 

1016 

1 123 


TAD 

LDIF 

1017 

3014 


DCA 

I DI F 

1020 

3055 


DCA 

55 

1021 

4531 


CALL MESAGE 

1022 

4543 

TEXT /%# 



1023 

4040 




1024 

4040 




1025 

4040 




1026 

4024 

T 



1027 

4040 




1030 

4040 




1031 

4040 




1032 

4040 




1033 

4040 




1034 

4040 




1035 

4040 




1036 

4031 

Y 



1037 

4040 




1040 

4040 




1041 

4040 




1042 

4040 




1043 

4040 




1044 

4040 




1045 

3103 

YC 



1046 

0114 

AL 



1047 

0340 

C 



1050 

4040 




1051 

4040 




1052 

4040 




1053 

4040 




1054 

4040 




1055 

041 1 

DI 



1056 

0640 

F 




KINETICS PROGRAM 

/ADJUST INFINITY VALUE? 
/YES 

/NO* WRITE RESULTS 

/INDEX 
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1057 

4543 

%# 


1060 

0000 

/ 


1061 

4407 

L00P4, 

FLPT 

1062 

541 1 


FGET I IT 

1063 

0014 


OUTPUT 

1 064 

0000 


FEXT 

1 06 6 

453 1 


CALL MESAG;i 

1066 

4040 


TEXT / 

1067 

4000 

/ 


1070 

4407 


FLPT 

107 1 

5412 


FGET I IY 

1072 

0014 


OUTPUT 

1073 

0000 


FEXT 

1074 

4531 


CALL MESAGE 

1075 

4040 


TEXT / 

1076 

4000 

/ 


1077 

4407 


FLPT 

1 100 

5413 


FGET I ICALC 

1 101 

0014 


OUTPUT 

1 102 

0000 


FEXT 

1 103 

4531 


CALL MESAGE 

1 104 

4040 


TEXT / 

1 105 

4000 

/ 


1 106 

4407 


FLPT 

1 107 

5414 


FGET I IDIF 

1 1 10 

0014 


OUTPUT 

till 

0000 


FEXT 

1 1 12 

4531 


CALL MESAGE 

1113 

4045 


TEXT / % 

1114 

4300 

#/ 


1115 

201 1 


ISZ IT 

1 1 16 

201 1 


ISZ IT 

1 1 17 

201 1 


ISZ IT 

1 120 

2012 


ISZ IY 

1 121 

2012 


ISZ IY 

1 122 

2012 


ISZ IY 

1 123 

2013 


ISZ ICALC 

1 124 

2013 


ISZ ICALC 

1 125 

2013 


ISZ ICALC 

1 126 

2014 


ISZ IDIF 

1 127 

2014 


ISZ IDIF 

1 130 

2014 


ISZ IDIF 

1 131 

2142 


ISZ NEGN 

1 132 

5261 

/PLOT 

JMP L00P4 
RESULTS ON CRT 

1 133 

4531 


CALL MESAGE 

1 134 

4543 

TEXT /%# 

1 135 

1 106 

IF 


1 136 

4020 

P 
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1 137 

1417 

LO 

1 140 

2440 

T 

1 141 

1 123 

IS 

1 142 

4004 

D 

1 143 

0523 

ES 

1 144 

1 122 

IR 

1 145 

0504 

ED 

1 146 

5440 

* 

1 147 

2201 

RA 

1 150 

1 123 

IS 

1 1 51 

0540 

E 

1 152 

130 5 

KE 

1 1 53 

3140 

Y 

1 1 54 

6154 

l. 

1 155 

4020 

P 

1 156 

2205 

RE 

1 157 

2 32 3 

SS 

1 160 

4003 

C 

1161 

1716 

ON 

1 162 

241 1 

TI 

1 163 

1625 

NU 

1 164 

0540 

E 

1 165 

4543 


1 166 

0000 

/ 

1 167 

7402 


1 170 

4540 


1 17 1 

4537 


1 172 

5773 


1 173 

1200 

PG 


HLT 

CALL MINMAX 
CALL PLOT 
JMP I PG5 
1200 
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/PAGE 5 

0 F r 11\ 3 1 0 !\ 0 Ek K I N ET ICS 

PROGRAM 

1 177 

1250 

* 1200 



1 200 

4531 

I TIE ST, 

CALL MESAGE 


1201 

4543 

TEXT /%} 

J 


1202 

1 124 

I T 

* 


1203 

0 522 

ER 



1204 

40 1 1 

I 



1205 

2340 

S 



1206 

2410 

TH 



1207 

0540 

E 



1210 

1625 

NU 



1211 

1 502 

MB 



1212 

0 522 

ER 



1213 

4017 

0 



1214 

0640 

F 



1215 

1 124 

IT 



1216 

0 522 

ER 



1217 

0 124 

AT 



1220 

1 1 17 

10 



1221 

1623 

NS 



1222 

5440 

S 



1223 

1 124 

IT 



1224 

0 522 

ER 



1225 

407 5 

= 



1226 

4000 

/ 



1227 

4533 


CALL SICONV 


1230 

3141 


DCA ITER 


1231 

1 121 


TAD LY 


1232 

4535 


CALL FLOAT 


1233 

4407 


FLPT 


1234 

6072 


FPUT TEST 


1235 

0000 


FEXT 


1236 

7200 


CL A 


12 37 

3055 


DCA 55 

/RESET FL. PT. OUTPUT 

1240 

1 141 


TAD ITER 


1241 

7 500 


SMA 

/IS ITER NEGATIVE? 

1242 

5247 


JMP ITRATE 

/NO: CONTINUE ITERATH 

1243 

7402 


HLT 

/YES: HALT 

1244 

7200 


CLA 


1245 

3055 


DCA 55 


1246 

5774 


JMP I RESTART 

/RESTART PROGRAM 

1247 

4775 

ITRATE, 

CALL CHECK 


1250 

7200 


CLA 

/GET DEV IN OBS QUANT 

1251 

1 122 


TAD LCALC 

/AND RMS ERROR 

1252 

3013 


DCA I CALC 


1253 

1 123 


TAD LDIF 


1254 

3014 


DCA I DIF 


1255 

1 124 


TAD LABS 


1256 

3015 


DCA IABS 
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1257 

1117 

TAD 

N 

1260 

7041 

CIA 


1261 

3142 

DCA 

NEGN 

1262 

4407 

FLPT 


1*26 3 

5150 

FGET 

ZERO 

1264 

6111 

FPUT 

ERROR 

1265 

0000 

FEXT 


1266 

4407 L00P5, FLPT 


1267 

5413 

FGET 

I I CALC 

1270 

0006 

FEXP 


127 1 

3075 

FMPY 

FACTOR 

1272 

1 126 

FADD 

FIN 

127 3 

6413 

FPUT 

I I CALC 

1274 

5415 

FGET 

I IABS 

1275 

2413 

F SUB 

I I CALC 

1276 

6414 

FPUT 

I IDIF 

1277 

3414 

FMPY 

I IDIF 

1 300 

1 1 1 1 

FAD D 

ERROR 

1301 

6111 

FPUT 

ERROR 

1 302 

0000 

FEXT 


1303 

2013 

I SZ 

I CALC 

1304 

2013 

I SZ 

I CALC 

1305 

2013 

I SZ 

I CALC 

1 306 

2014 

I SZ 

IDIF 

1307 

2014 

I SZ 

IDIF 

1310 

2014 

ISZ 

IDIF 

131 1 

2015 

ISZ 

I ABS 

1 312 

2015 

I SZ 

I ABS 

1313 

2015 

ISZ 

I ABS 

1314 

2142 

I SZ 

NEGN 

1315 

5266 

JMP 

LOOP5 

1316 

4531 

CALL 

MESAGE 

1317 

4543 TEXT 

/%# 


1 320 

2410 TH 



1321 

0540 E 



1322 

2215 RM 



1323 

2340 S 



1324 

0522 ER 



1325 

2217 RO 



1326 

2240 R 



1327 

0201 BA 



1330 

2305 SE 



1331 

0440 D 



1 332 

1716 ON 



1333 

4017 0 



1334 

0223 BS 



1335 

0522 ER 



1336 

2605 VE 



1337 

0440 D 




/INDEX 
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1 340 

2125 

QU 


1 341 

0116 

AN 


1342 

2456 

T. 


1343 

401 1 

I 


1 344 

2340 

5 


1 345 

0000 

/ 


1 346 

4407 


FLPT 

1 347 

5111 


FGET ERROR 

1350 

4114 


FDIV G 

1 351 

0002 


FSQRT 

1352 

6111 


FPUT ERROR 

1 353 

0014 


OUTPUT 

1 354 

0000 


FEXT 

1 355 

4531 


CALL MESAGE 

1 356 

4543 


TEXT /%# 

1357 

2410 

TH 


1 360 

0540 

E 


1361 

1 1 16 

IN 

. 

1 362 

0611 

FI 


1 363 

161 1 

NI 


1364 

2431 

TY 


1 365 

4026 

V 


1366 

0114 

AL 


1367 

2505 

UE 


1370 

4011 

I 


1371 

2340 

S 


1 372 

0000 

/ 


137 3 

5376 


JMP .+3 

1374 

0200 

RESTART, 

200 

1 37 5 

1637 

CHECK, 

CHK 

1 376 

4407 


FLPT 

1 377 

5126 


FGET FIN 
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/PG 6 OF FIRST ORDER KINETICS PROGRAM 
*1400 


1 400 

0014 


OUTPUT 

1401 

0000 


FEXT 

1 402 

7200 


CL A 

1 403 

1 147 


TAD ONE 

1 404 

7041 


CIA 

1 405 

1 141 


TAD ITER 

1 406 

3141 


DCA ITER 

1 407 

1141 


TAD ITER 

1410 

7440 


SZA 

1411 

577 1 


JMP I PG7 

1412 

4531 


CALL MESAGE 

1413 

4543 


TEXT /%# 

1414 

4040 



1415 

4040 



1416 

4024 

T 


1417 

1115 

IM 


1 420 

0 540 

E 


1 421 

4040 



1 422 

4040 



1423 

0000 

/ 


1 424 

7200 


CLA 

1 425 

1 125 


TAD M 

1 426 

1143 


TAD M2 

1427 

7550 


SPA SNA 

1 430 

5243 


JMP ALTERN 

1431 

4531 


CALL MESAGE 

1432 

4040 


TEXT / 

1 433 

4040 



1 434 

4040 



1 435 

2617 

VO 


1436 

1425 

LU 


1437 

150 5 

ME 


1 440 

4040 



1441 

4000 

/ 


1 442 

5254 


JMP LIST 

1443 

4531 

ALTERN, 

CALL MESAGE 

1444 

4040 


TEXT / 

1 445 

4040 



1 446 

0102 

AB 


1 447 

2317 

SO 


1 450 

2202 

RB 


1451 

01 16 

AN 


1452 

0305 

CE 


1453 

4000 

/ 


1 454 

4531 

LIST, 

CALL MESAGE 

1455 

4040 


TEXT / 

1 456 

4040 
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1 457 

4040 



1460 

4040 



1 46 1 

030 1 

CA 


1 462 

1403 

LC 


1463 

4040 



1 464 

4040 



1465 

4040 



1466 

4040 



1467 

4040 



1470 

4040 



147 1 

041 1 

DI 


1 472 

0640 

F 


1 473 

4543 

%# 


1474 

0000 

/ 


1475 

7200 


CLA 

1 476 

1 120 


TAD LT 

1 477 

301 1 


DCA IT 

1 500 

1 124 


TAD LABS 

1 501 

3015 


DCA IABS 

1 502 

1 122 


TAD LCALC 

1 503 

3013 


DCA I CALC 

1 504 

1 123 


TAD LDIF 

1 505 

3014 


DCA I DIF 

1 506 

1117 


TAD N 

1 507 

7041 


CIA 

1510 

3142 


DCA NEGN 

1 51 1 

3055 


DCA 55 

1512 

4407 

L00P6, 

FLPT 

1 513 

541 1 


FGET I IT 

1514 

0014 


OUTPUT 

1 515 

0000 


FEXT 

1516 

4531 


CALL MESAGE 

1517 

4040 


TEXT / 

1 520 

4000 

/ 


1 521 

4407 


FLPT 

1 522 

5415 


FGET I IABS 

1 523 

0014 


OUTPUT 

1 524 

0000 


FEXT 

1 525 

4531 


CALL MESAGE 

1 526 

4040 


TEXT / 

1 527 

4000 

/ 


1530 

4407 


FLPT 

1 531 

5413 


FGET I I CALC 

1 532 

0014 


OUTPUT 

1 533 

0000 


FEXT 

1 534 

4531 


CALL MESAGE 

1 535 

4040 


TEXT / 

1 536 

4000 

/ 


1 537 

4407 


FLPT 


/WRITE OUT RESULTS 
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1 540 

5414 


FGET 

I IDIF 

1 541 

0014 


OUTPUT 

1 542 

0000 


FEXT 


1 543 

4531 


CALL 

MESAGE 

1 544 

4045 


TEXT 

/ % 

1 545 

4300 

#/ 



1 546 

2011 


ISZ 

IT 

1 547 

201 1 


ISZ 

IT 

1 550 

201 1 


I SZ 

IT 

1 551 

2015 


ISZ 

I ABS 

1 552 

2015 


ISZ 

I ABS 

1 553 

2015 


ISZ 

I ABS 

1 554 

20 13 


I SZ 

I CALC 

1 555 

2013 


I SZ 

I CALC 

1 556 

2013 


I SZ 

I CALC 

1 557 

2014 


I SZ 

I OIF 

1 560 

2014 


I SZ 

IDIF 

1 561 

2014 


ISZ 

IDIF 

1 562 

2142 


ISZ 

NEGN 

1 563 

5312 


JMP 

L00P6 

1 564 

7402 


HLT 


1 565 

7200 


CL A 


1 566 

3055 


DCA 

55 

1 567 

5770 


JMP 

I RESTRT 

1 570 

1 133 

RESTRT* 

1133 

1571 

1600 

PG7 » 

1600 
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/PAGE 7 

OF FIxST OkDER 

KINETICS PROGRAM 



* 1600 



1600 

4407 

MODI NF.» 

FLPT 

/MODIFY INFINITY VALUE 

160 1 

5072 


FGET TEST 

/AND CONTINUE 

1602 

2111 


FSUB ERROR 


1603 

6 1 56 


FPU! TEM2 


1604 

0000 


FEXT 


1605 

7200 


CLA 


1606 

1 157 


TAD TEM2+1 


1607 

7510 


SPA 

/HAS RMS DEV DECREASED? 

1610 

5221 


JMP CHGSGN 

/NO: MOVE IN OPP. DIRECT 

16 11 

4407 


FLPT 

/YES: CONTINUE 

1612 

5111 


FGET ERROR 


1613 

6072 


FPUT TEST 


1614 

5126 


FGET FIN 


1615 

3067 


FMPY CORR 


1616 

6126 


FPUT FIN 


1617 

0000 


FEXT 


1620 

5705 


JMP I RUN 


1621 

7200 

CHGSGN* 

CLA 


1622 

4407 


FLPT 


1623 

5302 


FGET FL1 


1624 

4067 


FDIV CORR 


1625 

2302 


F SUB FL1 


1626 

3277 


FMPY HALF 


1627 

1302 


FADD FL1 


1630 

6067 


FPUT CORR 


1631 

3126 


FMPY FIN 


16 32 

6126 


FPUT FIN 


1633 

5111 


FGET ERROR 


1634 

6072 


FPUT TEST 


1635 

0000 


FEXT 


1636 

5705 


JMP I RUN 


1637 

0000 

CHK, 

0 


1640 

7450 


SNA 


1641 

527 5 


JMP AX 


1642 

4531 


CALL MESAGE 


1643 

4543 

TEXT /%£ 



1644 

1525 

MU 



1645 

1424 

LT 



1646 

1 120 

IP 



1647 

1431 

LY 



16 50 

4026 

V 



1651 

061 1 

FI 



16 52 

1640 

N 



1653 

0231 

BY 



1654 

4000 

/ 



1655 

4407 


FLPT 


1656 

0013 


INPUT 



29 





1657 

6067 


1 660 

0000 


1661 

4531 


1662 

4543 

TEXT 

1663 

261 1 

VI 

1664 

1606 

NF 

1665 

4075 

= 

1 666 

4000 

/ 

1667 

4407 


1670 

0013 


1671 

6 126 


1672 

0000 


1673 

2141 


1 674 

5637 


1675 

2141 

AX, 

1676 

5705 


1677 

7777 

HALF 

1 700 

3777 


1701 

7774 


1702 

0000 

FL1, 

1703 

3777 


1704 

7774 


1705 

0563 

RUN, 


FPUT CORK 
FEXT 

CALL MESAGE 

/%# 


FLPT 
INPUT 
FPUT FIN 
FEXT 

ISZ ITER 
RETURN CHK 
ISZ ITER 
JMP I RUN 
FLTG 0.5 


FLTG 1.0 


Y CALC 







ALTERN 

1 443 

AX 

1675 

A 1 

0622 

B 

0103 

CALL 

4400 

CHECK 

137 5 

CHGSGN 

1621 

CHK 

1637 

CONST 

0144 

CORK 

0067 

DECPRT 

0132 

ERROR 

0111 

FACTOR 

0075 

FEXP 

0006 

FIN 

0126 

FIX 

0134 

FLAG 

0153 

FLOAT 

0135 

FLOG 

0007 

FLPT 

4407 

FL1 

1702 

FSQRT 

0002 

G 

0 114 

HALF 

1677 

I ABS 

0015 

I CALC 

0013 

IDIF 

0014 

INPUT 

0013 

IT 

001 1 

I TER 

0141 

ITRATE 

1247 

ITTEST 

1200 

IY 

0012 

LABS 

0124 

LCALC 

0122 

LDIF 

0123 

LIST 

1454 

LOOP 1 

0532 

L00P3 

0633 

L00P4 

106 1 

LOOPS 

1266 

L00P6 

1512 

LSG 

0136 
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LT 

0120 

LY 

0121 

M 

0125 

MAX 

0172 

MESAGE 

0131 

MIN 

0167 

MINMAX 

0140 

MODINF 

1600 

M2 

0143 

N 

01 17 

NEG 

0627 

NEGN 

0142 

NEXT 

0430 

ONE 

0147 

OUTPUT 

0014 

PG2 

0377 

PG3 

0577 

PG4 

0772 

PG5 

1 173 

PG7 

1 57 1 

PLOT 

0137 

R 

0106 

READ 

0410 

RESTAR 

1374 

RESTRT 

1570 

RETURN 

5400 

RUN 

1705 

S 

0100 

SICONV 

0133 

START 

0200 

TEMP 

0175 

TEM2 

0156 

TEST 

0072 

XLOC 

0154 

XMIN 

0156 

XRNG 

0161 

XZERO 

0172 

Y CALC 

0563 

YLOC 

0155 

YMIN 

0164 

YRNG 

0167 

YZERO 

0175 

ZERO 

0150 















9 


4400 0000 

4401 7300 

4402 7404 

4403 0332 

4404 7440 

4405 5210 

4406 2200 

4407 5600 

4410 7201 

4411 3326 

4412 1333 

4413 3011 

4414 1334 

4415 3012 

4416 1117 

4417 7041 

4420 3327 

4421 2327 

4422 4407 

4423 5411 

4424 6172 

4425 5412 

4426 6175 

4427 0000 

4430 2011 

4431 2011 

4432 2011 

4433 2012 

4434 2012 

4435 2012 

4436 4407 

4437 5411 

4440 2172 

4441 6161 

4442 5412 

4443 2175 

4444 6167 

4445 0000 


/SUBROUTINE FOR OBTAINING MINIMUM AND MAXIMUM 
/VALUES OF ARRAYS 
/SUBROUTINE MINMAX 

/SPECIAL VERSION FOR FIRST ORDER KINETICS PROGRAM 

RETURN= 5400 

IX= 1 1 

I Y= 12 

N= 1 17 

TEMX=172 

TEMY=175 

TEMX 1= 161 

TEMY1=167 

XMIN=156 

XRNG= 161 

YMIN=164 

YRNG=167 

FLPT= 4407 

*4400 

MINMAX, 0 

CLA CLL 
OSR 

AND MASK 
SZA 

JMP .+3 
ISZ MINMAX 
RETURN MINMAX 
CLA I AC 
FLAG 
XLOC 
IX 

YLOC 
IY 
N 


START, 


LOOP/ 


DCA 
TAD 
DCA 
TAD 
DCA 
TAD 
CIA 
DCA 
I SZ 
FLPT 
FGET 
FPUT 
FGET 
FPUT 
FEXT 
ISZ IX 
I SZ 
I SZ 
I SZ 
I SZ 
I SZ 
FLPT 
FGET 
F SUB 
FPUT 
FGET 
F SUB 
FPUT 
FEXT 


/FLAG = +1 
/SET INDEX 


VALUES 


NEGN 

NEGN 

I IX 
TEMX 
I IY 
TEMY 


IX 

IX 

IY 

IY 

IY 

I IX 
TEMX 
TEMX 1 
I IY 
TEMY 
TEMY 1 


/FIRST VALUE OF X 
/SAVE IT 

/FIRST VALUE OF Y 
/SAVE IT 

/GET MINIMUM VALUES 
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4446 

4447 

4450 

4451 

4452 

4453 

4454 

4455 

4456 

4457 

4460 

4461 

4462 

4463 

4464 

4465 

4466 

4467 
4470 
447 1 

4472 

4473 

4474 

4475 

4476 

4477 
4 500 

4501 

4502 

4503 

4504 

4505 

4506 
4 507 

4510 

4511 

4512 

4513 

4514 

4515 

4516 

4517 
4520 
4 521 
4 522 

4523 

4524 

4525 

4526 
4 527 

4530 

4531 

4532 

4533 

4534 


7200 


CL A 

1 162 


TAD TEMX1 

7 500 

Al* 

SMA 

5256 


JMP NEXT 

4407 


FLPT 

541 1 


FGET I IX 

6172 


FPUT TEMX 

0000 


FEXT 

7200 

NEXT* 

CL A 

1 170 


TAD TEMY1 

7500 

A2* 

SMA 

5266 


JMP NEXT1 

4407 


FLPT 

5412 


FGET I IY 

6175 


FPUT TEMY 

0000 


FEXT 

2327 

NEXT1* 

ISZ NEGN 

5230 


JMP LOOP 

7240 

' 

CLA CM A 

1326 


TAD FLAG 

7440 


SZA 

531 1 


JMP MAX 

4407 


FLPT 

5172 


FGET TEMX 

6156 


FPUT XMIN 

5175 


FGET TEMY 

6164 


FPUT YMIN 

0000 


FEXT 

2326 


ISZ FLAG 

7200 


CLA 

1330 


TAD CHNG 

3250 


DCA A1 

1330 


TAD CHNG 

3260 


DCA A2 

5212 


JMP START 

7200 

MAX* 

CLA 

1331 


TAD NORM 

3250 


DCA A1 

1331 


TAD NORM 

3260 


DCA A2 

4407 


FLPT 

5172 


FGET TEMX 

6161 


FPUT XRNG 

5175 


FGET TEMY 

2164 


FSUB YMIN 

6167 


FPUT YRNG 

0000 


FEXT 

5600 


RETURN MIf 

0000 

FLAG* 

0 

0000 

NEGN* 

0 

7510 

CHNG, 

SPA 

7 500 

NORM, 

SMA 

2000 

MASK* 

2000 

3200 

XLOC* 

3200 

2400 

YLOC* 

2400 


/STORE MIN VALUES 


/NOW GET MAX VALUES 


/STORE RANGES 


/XMIN ASSUMED TO BE ZERO HERE 
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✓ * 


/PLOT DATA; XMIN IS ASSUMED TO BE ZERO 
/YRANGE IS ASSUMED TO INCLUDE ZERO 
*2000 


2000 

0000 

PLOT, 

0 

2001 

7200 


CL A 

2002 

1777 


TAD YLOC 

2003 

3012 


DCA IY 

2004 

1776 


TAD XLOC 

2005 

3011 


DCA IX 

2006 

1363 


TAD YFLOC 

2007 

3365 


DCA JY 

2010 

1362 


TAD XFLOC 

2011 

3364 


DCA JX 

2012 

1 117 


TAD N 

2013 

7041 


CIA 

2014 

3361 


DCA INDEX 

2015 

4407 

LOOPA, 

FLPT 

2016 

541 1 


FGET I IX 

2017 

2156 


FSUB XMIN 

2020 

4161 


FDIV XRNG 

2021 

3346 


FMPY CONST 

2022 

0000 


FEXT 

2023 

47 57 


CALL FIX 

2024 

3764 


DCA I JX 

2025 

4407 


FLPT 

2026 

5412 


FGET I IY 

2027 

2164 


FSUB YMIN 

2030 

4167 


FDIV YRNG 

2031 

3346 


FMPY CONST 

2032 

0000 


FEXT 

2033 

47 57 


CALL FIX 

2034 

3765 


DCA I JY 

2035 

2011 


ISZ IX 

2036 

201 1 


ISZ IX 

2037 

2011 


ISZ IX 

2040 

2012 


ISZ IY 

2041 

2012 


ISZ IY 

2042 

2012 


ISZ IY 

2043 

2364 


ISZ JX 

2044 

2365 


ISZ JY 

2045 

2361 


ISZ INDEX 

2046 

5215 


JMP LOOPA 

2047 

4407 


FLPT 

20 50 

5351 


FGET ZERO 

2051 

2164 


FSUB YMIN 

2052 

4167 


FDIV YRNG 

2053 

3346 


FMPY CONST 

2054 

0000 


FEXT 

2055 

47 57 


CALL FIX 

2056 

3366 


DCA YY 


/FIX POINTS 


/GET AXIS LOCATIONS 
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20 57 

3355 


DCA 

X VAL 

2060 

7200 

STARTS, 

CLA 


2061 

1366 


TAD 

YY 

2062 

3356 


DCA 

Y VAL 

2063 

1354 


TAD 

MAX V 

2064 

3360 


DCA 

I ND2 

2065 

1356 


TAD 

Y VAL 

2066 

6063 


DYL 


2067 

1355 


TAD 

X VAL 

2070 

6057 


DXS 


207 1 

2355 


ISZ 

X VAL 

2072 

2355 


ISZ 

X VAL 

2073 

2360 


ISZ 

I ND2 

2074 

5267 


JMP 

.-5 

2075 

7200 


CLA 


2076 

1354 


TAD 

M AX V 

2077 

3360 


DCA 

I ND2 

2 100 

3355 


DCA 

X VAL 

2 101 

3356 


DCA 

YVAL 

2102 

6053 


DXL 


2103 

1356 


TAD 

Y VAL 

2104 

6067 


DYS 


2105 

2356 


ISZ 

YVAL 

2106 

2356 


I SZ 

YVAL 

2 107 

2360 


ISZ 

IND2 

2110 

5303 


JMP 

.-5 

2 111 

7200 


CLA 


2 1 12 

1363 


TAD 

YFLOC 

2 113 

3365 


DCA 

JY 

21 14 

1362 


TAD 

XFLOC 

2115 

3364 


DCA 

JX 

2116 

1117 


TAD 

N 

2117 

7041 


CIA 


2120 

3360 


DCA 

I ND2 

2 121 

7200 

LOOPB, 

CLA 


2 122 

1764 


TAD 

I JX 

2123 

6053 


DXL 


2124 

7200 


CLA 


2125 

1765 


TAD 

I JY 

2126 

6067 


DYS 


2127 

7200 


CLA 


2130 

1367 


TAD 

M 10 

2131 

3361 


DCA 

INDEX 

2132 

2361 


I SZ 

INDEX 

2133 

5332 


JMP 

1 

2134 

2364 


I SZ 

JX 

2 135 

2365 


ISZ 

JY 

2 136 

2360 


I SZ 

IND2 

2 137 

5321 


JMP 

LOOPB 

2140 

7200 


CLA 



/DRAW AXES 


/X AXIS 


/Y AXIS 


/AXES DRAWN 


/PLOT POINTS 


/POINT PLOTTED 


/DELAY 
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